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Abstract 

We study the moment space corresponding to matrix measures on the unit circle. Mo- 
ment points are characterized by non-negative definiteness of block Toeplitz matrices. This 
characterization is used to derive an explicit representation of orthogonal polynomials with 
respect to matrix measures on the unit circle and to present a geometric definition of canon- 
ical moments. It is demonstrated that these geometrically defined quantities coincide with 
the Verblunsky coefficients, which appear in the Szego recursions for the matrix orthogonal 
polynomials. Finally, we provide an alternative proof of the Geronimus relations which is 
based on a simple relation between canonical moments of matrix measures on the interval 
[-1,1] and the Verblunsky coefficients corresponding to matrix measures on the unit circle. 
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1 Introduction 

In recent years considerable interest has been shown in moment problems, orthogonal polynomials, 
continued fractions and quadrature formulas corresponding to matrix measures on the real line 
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or on the unit circle. Early work dates back to iKreinI (11949), while more re c ent resu l ts on matrix 



mea sures on the real line can be found in the papers of 
and 



Defez et al. 



Rodman! fll990h . lOuranl fll995l . Il996f ) 



(I2OOOI ) among many others. Additi o nally, several aut h ors have discussed matrix 



measures on the unit circle \see 



Delsarte et al. 



(|l978h.lGeronimd (|l98ll).lMarcellan and Rodriguez 



1989), ISinap and Van Asschel (11994 Il99fil ) 



Yakhlef and Marcellan 



2001 



20021 ). 



Cantero et al. 



(12003h ]. 



The purpose of the present paper is to investigate some geometric properties of the moment space 
corresponding to matrix measures on the unit circle. In Section 2 we present a characterization 
of the moment space in terms of nonnegative definiteness of block Toeplitz matrices. We also 
provide a geometric definition of canor iical moments of matrix rn easures on the unit circle, which 
generalizes the scalar case discussed by lDette and StuddenI (119971 ) in a nontrivial way. In Section 3 
an explicit determinantal representation of orthogonal matrix polynomials with respect to matrix 
measures on the unit circle is presente d , whic h generalizes the classical representation in the one- 
dimensional case [see e.g. iGeronimusI (1l962l )]. These results are used to identify the canonical 
moments as Verblunsky coefficients, which appear in the Szego relation s for the corresponding^ 
orthon or mal and reversed matrix polynomials [see 



Delsarte et al. 



(1l978l ). ISinap and Van Assche 



(119961 ) or iDamanik et al.l (j2008l )]. In particular our results provide a geometric definition of 
Verblunsky coefficients corresponding to matrix measures on the unit circle. Roughly speaking, 
the Verblunsky coefficient of order m can be characterized as the distance of the mth trigonometric 
moment to a center of a matrix disc relative to the diameter of this disc (see Section 3 for 
more details). Finally, in Section 4 these results are used to present an alternative proof of the 
Geronimus relations for monic orthogonal polynomials, which describe the relation between the 
coefficients in the three-term recursive relation of orthogonal polynomials with respect to a matrix 
measure on a compact interval and the coefficients in the Szego recursion of an associated matrix 
measure on the unit circle. 



2 The moment space of matrix measure on the unit circle 

A matrix measure fi on the unit circle is defined as apxp matrix of (real valued) Borel measures 
/i = (/iij)i,j=i,...,p on the unit circle d3 = {2 G C| \z\ = 1} such that for each Borel set A C d3 
the matrix fi{A) is nonnegative definite, i.e. /^(A) > 0. Throughout this paper we use the usual 
parametrization z = e*^, 6 G [— vr, vr) and the notation fi{6) for the sake of simplicity. The kth 



2 



moment of a matrix measure fi on the unit circle is defined by 

/TT 

where = = f^^ cos {k9)dfi{9), (3^ = Pk{f^) fl^^^^ {k9)dfi{6) {k = 0,1,...) are the 

trigonometric moments and the dependence on the given measure fi is omitted in the no- 
tation, whenever it is clear from the context. Throughout this paper let m G No A(yu) = 
{ao,ai, Pi, . . . ,am, Pm) € (MP><P)2m+i (jgnote the vector of trigonometric moments of order m 
and define 

(2.2) M2m+i = {A(/x) I /i is a matrix measure on dB} C (Mpxp)2"^+i 

as the (2m + l)th moment space of matrix measures on the unit circle. The set M.2m+i and its 
interior Int(A^2m+i) can be characterized as follows. 

Theorem 2.1 A = (ao, ai, . . . , a^, Prn) e M2m+i if and only if 

m m 

(2.3) Yl trace{B,B*T,.,) > V Bq, . . . , 5„ G C^^^, 

i=0 j=0 

where the matrices r_m, T^m+i, ■ ■ ■ , Tm are defined in h2.1\) . 

A = (ao, ai, . . . , am, Pm) G Int{M.2m+i) if and only if there is strict inequality in h2.3^) except 
zfBo = --- = Bm = 0. 



Proof: We start with a proof of the first part. Assume that A G M.2m+i and consider matrices 
Bo,...,Bme CP^P. With the notation 

m 

(2.4) B{e) = Y,Bke''' ^ee[-n,7r)) 

k=0 

it follows that the polynomial P{6) = B{6){B{6))* is obviously nonnegative definite, i.e. 

(2.5) P{9) = B{e){B{e))* > G [-TT, tt). 

A straightforward calculation shows that the polynomial P can be represented as 

m 

(2.6) P{e) =Do + J2^k cos {ke) + Ek sin {kO), 

k=l 
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where the p x p matrices Dq, ■ ■ • , -D^, -Ei, • • • , are defined by Dq = Aq, and for k = 1, . . . ,m 



and 



m—k 



Ak = J2 ^k+iB! and A_k = A^ 



1=0 



Because it is easy to see that the moment space A42m+i is the convex hull of the set 

l^aa*, cos (6')aa*, sin (6')aa*, cos (m6')aa*, sin (m6')aa*) a G C^, 9 E [— 7r,7r)|, 



a similar argument as in Corollary 2.2 of iDette and Studded (120021 ) now shows that (12.51) and 
dMD imply 



< trace(Do«o) + trace (Dfea/c) + tTa,ce{EkPk 

k=l 

= trace] / d{Dofi{e)) + ^ cos {k9)d{Dkfxi9)) + j sm{ke)d{Ekfi{e)) 

m ^ 

trace I / ^ e'''^d{Akfi{e)) 

(njY lib I III, IX. \ nq^ lib I III, Ix. 

•'-^k=o \«=o / ■^-^k=\ \/=0 , 

EE / e'^''''^'d{BkB*fi{e)) 
k=o 1=0 ) 

m m 
k=0 1=0 

which proves (12.31) . On the other hand assume that the inequality (12. 3p is satisfied for all matrices 
Bq, . . . , Bm G C^^P and consider a nonnegative definite matrix polynomial 



m—k 



(2.7) 



P{e) = Do + ^Dk COS {ke) + Ek sin (kO) > G [-vr, vr) 



fc=i 



with hermitian matrices Do, . . . , Ei, . . . , Em G C^^^. It now follows from iMalyshevI (Il982l ) 
that there exists a matrix polynomial 



B{d) = J2Bke''' 



k=0 

4 



such that P{0) = B{6){B{6))* , and the same calculation as in the first part of the proof yields 

m mm 

trace(DoCto) + trace(Dfcttfc) + trace{EkPk) = trace(-Bii?*rj_j) > 0. 

i=0 j=0 



k=l 



By similar arguments as in Lemma 2.3 of iDette and Studded (120021 ) it follows that this is suffi- 
cient for A G M.2m+1 ■ 

Finally, the second part of the Theorem is shown similarly observing the fact that 
(ao, tti, . . . , a™, (3jn) e Int{M2m+i) if and only if 

m 

trace(Dotto) + trace(-Dfcafc) + tTace{EkPk) > 

k=l 

for any nonnegative definite polynomial P{0) of the form (12.61) with P (6) ^ G [— tt, tt ) . This 
characterization can be shown by the same arguments as presented in iDette and Studderu (120021 ) 
who proved a corresponding statement for the moment space of matrix measures on the interval 
[0,1]. □. 



Throughout this paper let 

(2.8) = r„(/i) 



r \ 



To J 



m+l) xp(m+l) 



denote the Block Toeplitz matrix, where the blocks Fj = Fi(yu) {i = — m, . . . , m) are the moments 
of a matrix measure fi on the unit circle defined by (12. ip (note that is hermitian). The 
following characterization of the moment space Al2m+i by nonnegative definiteness of Toeplitz 
matrices is now easily obtained. 



Corollary 2.2 Assume that X = (ao, ai. Pi, ... , am, Pm] 
Ii2.8\) with Tk = ak + iPk and F_fc = - iPk- 

(a) \ E M.2m+i if and only ifT^ > 0. 

(b) X E Int{M.2m+i) if and only if Tm > 0. 



e (Mpxp)2™+i and that is defined by 



Proof: We only proof part (a); part (b) is shown by similar arguments. First assume that 
XeM 2m+i) then we obtain from Theorem 12JJ for all matrices Bq, . . . , B.^ G C'p^^ 



m m 



2^2^tTace{B,B*Tj_,) > 0. 

8=0 j=0 
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Consequently, if ao, . . . , G C^, a = (a^, . . . , a^)^ G Cp(™+i) we put Bi = (a^, 0, . . . , 0) G C^^p 
= 0, . . . , m) and it follows 

mm mm 

a*Tma = trace(aa*Tm) = trace(aia*rj_j) 

1=0 j=0 1=0 j=0 

which shows that the matrix is nonnegative definite. To prove the converse assume that 
Tm > 0, i.e. 



:2.9) 



< a*TmCi = trace(aja*r 

1=0 j=0 



for all a = (a^, . . . , a^f G C^^^+i). If 5o, . . . , 5„ G C^^^, and aj'^ denotes the ith column of the 
matrix Bj {j = 0, . . . ,m, i = 1, . . . ,p), then 



1=1 

and we obtain from (12.91) 

mm p m m 

5^^trace(S,S;r,_,) = j^J^J^ trace (af^ (af )*r,_i) > 0. 



1=0 j=0 



k=l i=0 j=0 



By Theorem 12.11 it follows that A G M.2m+i, which completes the proof of the Corollary. □ 



With the aid of Theorem 12.11 and Corollary 12.21 we are now able to define geometrically canonical 
moments for matrix measures on the unit circle. It turns out that these geometrically defined 
quantities are exactly the Verb 
duced by 



Damanik et al. 



unsky coefficients of matrix measures on the unit circle as intro- 
(120081 ) (see Section 3 where we prove this identity). For this purpose 



let W denote a p x p matrix and define 



:2.io) 



A = A{W) 



\ 



F F 

— m — m+l 

w* r_™ 



Tm W \ 

Tm— 1 Tm 

To Ti 

r-i To / 



m+2)xp(m+2) 



Let r*^™) = (r_m, r_m+i, . . . , r^-i, r^) G (CP^P)^™+^ denote a vector of moments of a matrix 
measure on the unit circle, that is (oq, ai, Pi, . . . , am, Pm) ^ M.2m+ii where = + iPk- Define 
Pp{m) as the set of all matrix measures on the unit circle with moments of order m given by F*^™\ 
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that is Tj = f^^ e"^^ djjiiO) (j = —m, . . . , m). By Corollary 12.21 it follows that the matrix W is the 
(m + 2)th moment of a matrix measure fi G Ppf™) if ^^"^ only if yl(iy) > 0. We assume without 
loss of generality that (ao,ai,/3i, . . . , am., Pm) £ IntfA^ 2^ 4 -1 ) wh ich is equivalent to > by 
Corollary 12.21 . From Theorem 1 in iFritzsche and KirsteinI (119871 ) it follows that 



A{W) > 

if and only if there exists a p x p matrix U with UU* < Ip such that the matrix W can be 
represented as 



f2.111 



W 



(Fi . . . Fm) T^-i (r_m • • • r_i)* + Ll/^'^u Rll^, 



where the matrices and Rm are defined by 



(2.12) 
(2.13) 



Rr, 



ro — (ri . . . Fm) (Fi . . . F^)* , 
ro — (F_m • • • F_i) T^ii (F_m . . . F„i) 



respectively. Note that the matrices Lm and Rm are Schur co mplements of the po s itive definite 
matrix and as a consequence are also positive definite [see lHorn and Johnsohru (119851 )]. This 
means that that the matrix W is the (m + 2)th moment of the matrix measure ^ G V-p(ni), if and 
only if it is an element of the "ball" 

(2.14) Km CP><^|L-V2(iy _ Mm)R~^"^ = U, UU* < /J , 
where the "center" of the ball is given by the matrix 

(2.15) M„ = (Fi . . . F^) T-ii (F_„ . . . F_i)* . 

We are now in a position to define the canonical moments of a matrix measure on the unit circle 
(or Verblunsky coefficients as shown in Section 3). 



Definition 2.3 Let jj, denote a matrix measure on the unit circle with moments F^ = + ijSk 
(k > 0), X2m+M = (ao,ai,/5i, • • .,am,f3m) G (Mpxp)-+i (m > 0) and define 

(2.16) N{fi) = min {m e N I \2m+M e dM2m+i} , 

as the minimum number m G N such that X2m+i is a boundary point of the moment space M.2m+i 
(if A2m+i £ Int{M.2m+i) for all m & N we put N{fi) = 00). For each m = 0, . . . , N{fi) — 1 the 
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quantity 

(2.17) Am+i=A^+^{^) = L;,'/^T^+^ - Mm) R~'^^ 

= [ro — (Fi, . . . , Fm) T^ii (Fi, . . . , r^)*] ^ 

X (Fm+1 — (Fl, . . . , Tm) T^-l (F_m, • • • , F_i) ) 

X [Fo - (F_^, . . . , F_i) T-ii (F_^, . . . , F_i)*] 
is called the {m + l)th canonical moment of the matrix measure fi. 



Definition 12.3) is a generalization o : 



the definition of canonical moments of scalar measures on the 



unit circle in iDette and StuddenI (jl997l ). In general the explicit representation of the canonical 
moments in terms of the moments Fq, Fi, ... is very difficult. For example if m = we have 

(2.18) A^ = Fo '/'FiFq 
and in the case m = 1 we obtain from Definition 12.31 

(2.19) A2 = (Fo - FiFo 'F_i)-'/' (F2 - FiFq ^Fi) (Fq - F^iFoFi)"^/' 

In the following section we will demonstrate that the quantities defined by Definition 12.31 are the 
well known Verblunsky coefficients, which are usually obtained from the recursive relations of 
the orthonorrnal poly nomials with respect to matrix measures on the unit circle [see for example 



Delsarte et al. 



(119781 ) where these matrices do not hav e any special na r ne. ISinap and Van Assche 



(Il996l ) where they are called reflection coefficients or iDamanik et al.l (120081 )]. For this purpose 
we use an explicit determinant representation of the matrix orthogonal polynomials, which is of 
own interest and given in the following section. 



3 Orthogonal matrix polynomials 

A p X p matrix polynomial is a p x p matrix with polynomial entries. It is of degree n if all the 
polynomial entries are of degree less than or equal to n and is usually written in the form 



(3.1) 



i=0 



with coefficients Ai G C^^^ and z G C. Recall that for matrix polynomials P and Q the right 
and left inner product are defined by 

Pie''ydMQ{e''), 



(3.2) 
(3.3) 



{P,Q)r 
{P,Q)l 



Pie'^)dfx{e)Q{e'' 
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respectively [see for example ISinap and Van Asschd (Il996l )]. The matrix polynomials P and Q 
are called orthogonal with respect to the right inner product (■, ■)r if 



(3.4) 



R 







and orthogonality with respect to the left inner product (■, ■)l is defined analogously. The matrix 
polynomials Pq{z), Pi{z), P2{z), . . . are called orthonormal with respect to the right inner product 
if for each m G Nq Pm{z) is of degree m, Pm{z) and P^'i^) ^-^e orthogonal with respect to (■, ■)r 
whenever m ^ m and 



(3.5) 



^Pmi Pm) R 



•-pi 



where Ip denotes the p x p identity matrix. Orthonormal polynomials with respect to the left 
inner product (■, ■)l are defined analogously. Orthonormal polynomials with respect to the inner 
products (■, ■)r and (■, ■)l are determined uniquely up to multiplication by unitary matrices. In 
the following discussion we will derive an explicit representation of these polynomials in terms of 
the moments of matrix measure which generalizes the well known determinant representation 
in the scalar case [see for example iGeronimud (119461 )]. 

For this purpose consider a matrix measure /i on the unit circle with moments r_m, . . . , and 
recall the definition of the corresponding block Toeplitz matrix in (12.81) . We define for m G N 
matrix polynomials by 



(3.6) 
(3.7) 



('^i(^))ij=i,..,p' 

(-^ii(^))ij=l,...,p ' 



where the elements T^{z) and TMz) in these matrices are given by the determinants 



(3.8) 



and 



(3.9) 



To 

r-i 

TlJz) 



ri,{z) 



To 



-m+2 



m—1 



To 

— m+l 



,P 



m—1 
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respectively, and the matrices Tj^^^^ (and f j^^^^,) are obtained replacing the jth row (and the 
zth column) in the matrix T^m+k by efz^ (and ejz"^~^). The following result shows that these 
polynomials are orthogonal with respect to the given matrix measure /z. 

Theorem 3.1 For a given matrix measure fi on the unit circle let "^^{z) and ^^^(-z) (m E N) 
denote the matrix polynomials defined by liS. 6\) and Ij^S. 7\ ), respectively, then we have 

(3.10) {z'lp,K)R = 0, {k = 0,...,m-l) ; vl/«)^ = |T^|/, 

{^t^'lp)L = 0, (A; = 0,...,m-1) ; {^H^, z"^I,)l = \T^\Ip. 

Proof: We will only give a proof for the polynomials \E'^(z), the remaining part of Theorem 13. II 
is shown similarly. The element B[j in the position {i,j) of the matrix 

B^:={z'l,^^)n= re-^''^M^)(^?(e^')),,=i,...,, (fc = 0,...,m), 

is given by 

(3.11) ^5 = E f e-^'^'Tl^{e^')dUe). 

An expansion of the determinant T^[e^^) with respect to the {mp + j)th row yields 

m 

^3 Ty[e^^) = ( — 1 ) ('^+")p+-?'+'e^"^ |2^("^p+i).("p+0| 

n=0 

where the matrix jg obtained from by deleting the (mp + j)th row and {np + l)th. 

column. If '^n,ij = Jlj^^^^^dfiij denotes the element of the matrix r„ in the position {i,j), where 
n G {—m, . . . , m}, it follows that 

m p 

(3.13) Bfj = ^(-l)('^+")P+^'+' |T^"'^+''^'^"''+'Vn-fc,*/. 

n=0 1=1 

Now it is easy to see that the right hand side of (13.131) is the determinant of the matrix Tm, where 
the {mp + j)th row has been replaced by the vector 

(T— • • • ) "^—kjipi T— • • • ■> T— fc+l,jp) • • • ■> ^m—l—k^ili • • • i ^m—l—k,ip^ Tm— fe,il; • • • ; '^m—k^ip) 

Consequently, if A; G {0, . . . , m — 1} the {mp + j)th and {kp + i)th row in this matrix coincide 
and we have B^^ = 0, which proves the first identity in (13.101) . 
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For a proof of the second identity we note that in the case k = m and i ^ j the same argument 
yields Bij = 0. \i k = m and z = j it follows that Bij is exactly the determinant of the matrix 
Tmi which completes the proof of the first assertion of Theorem 13.11 □ 

In the following discussion we derive several consequences of the representations (13.61) and fl3.7p . 
which will be useful to identify the canonical moments as Verblunsky coefficients. In particular 
we determine the corresponding leading coefficients and identify the orthonormal polynomials 
with respect to the measure ^. For this purpose recall that a matrix polynomial of the form (13.11) 
is called monic, if the coefficient of the leading term is the identity matrix, that is An = Ip. 

Corollary 3.2 For a given matrix measure fi on the unit circle let \l'm(^) o^nd "^mi^) defined 
by 1^3. 6i) and and consider for m < N{^) the matrix polynomials 



(3.14) $«(z) = ^^{z)\T„,\-'R^, 

(3.15) $^^(z) = \Trr,\-'L^^i{z), 

where the matrices Rm and Lm are defined by ( 12.13^ and f l2.12p . respectively. The polynomials 
^m(^) (and $^(z)j are monic orthogonal matrix polynomials with respect to the right (and left) 
inner product {■,-)r (and 
Similarly, define for m < N{fi) 

(3.16) 0^(z) = Ki^)\T^\-'Rll^ 

(3.17) 0^(z) = |r„|-iLV2vl>^(^), 

then the matrix polynomial (and (f)^{z)) are orthonormal polynomials with respect to the 

right (and left) inner product (■, (and (■, ■)l). The leading coefficients of (f)^{z) and 0^(-2) are 
given by Rm^'^ and Lm^"^ , respectively. 

Proof: In the first part we will prove that the leading coefficients of the polynomials '^^{z) and 
\I'^(z) defined by (13. 6p and (13. 7p are given by 

(3.18) = \Tm\R'^, 

(3.19) = |Tm|I^~\ 

respectively. With these representations we obtain from Theorem 13.11 

/ayfi vl>^\o = IT \{L^Y ■ /v|/^ \l/^\r = IT \iL^Y 

\^mi^m/R \-^m\\-^rnJ ) \^ mi ^ ml L \-^m\\-^m) i 
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and the assertion of the Corollary follows by a straightforward calculation. 

In order to prove (13.181) and fl3.19p we restrict ourselves to the first case; the second case is shown 
similarly. Observing the definition of the determinants T^{z) in fl3.8p we obtain for the entry in 
the position j) of the leading coefficient of the matrix polynomial '^^{z) 

where we have used an expansion of the determinant with respect to the {mp + j) row and the 
matrix jg obtained from T.^ by deleting the (mp + j)th row and (mp + z)th column. 

This means that [L^^. is the entry in the position {mp + i, mp + j) of the adjoint of the matrix 
Tm {hi = 1, • • • and consequently L^/|Tm| is the p x p block in the position (m + 1, m + 1) 
of the matrix T~^, which is given by 



(ro-(r_^...r_i)r„ii(r_™...r_i)*) 



-1 



[see e.g. 



Horn and JohnsohnI (119851 )]. This proves the assertion (13.181) and completes the proof 
of the Corollary. □. 



We are now in a position to identify the canonical moments introduced in Definition 12.31 as 
Verblunsky coefficients which are defined as coefficients in the Szego relation of the matrix or- 
thonormal polynomials 0^(-2) and (pni^)- this purpose we introduce for a given matrix 

polynomial P„ of degree n the corresponding reversed polynomial 

r ■ 

where z denotes the complex conjugation of z G C. Obviously we have for any p x p matrix A 



PJz) 



/IP„(^) = Pn{z)A*. 



In the following discussion let k 



R 



Rm ^ and 



Lm^'^ {m 



, N{iJi) — 1) denote the 



leading coefficients of the orthonormal matrix polynomials (f)i^{z) and (f)j^{z) with respect to the 



b^{z) and 0^1 

right and left inner product induced by the matrix measure fi and define the matrices 
(3.20) P™=K+i)"''t^ and = {m = 1, . . . , N{fi) - 1). 

Then it follows from 



Damanik et al 



(120081 ) that there exist p x p matrices such that the 
orthonormal matrix polynomial with respect to the measure /i on the unit circle satisfy the Szego 
recursions 



(3.21) 
(3.22) 



-Pm0m+l(2;) = Hm+l(p^{z) 



'm+1 • 
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The matrices are uniquely determined and called Verblunsky o r reflection coefficients, b ecause 
they were introduced for the scalar case in two seminal papers by I VerblunskyI (jl935l . Il936l ). The 
flnal result of this section shows that the Verblunsky coefficients coincide with the canonical 
moments introduced in Deflnition 12.31 



Theorem 3.3 Let fi denote a matrix measure on the unit circle and assume that < m < N{fi). 
IfAm+i is the {m+l)th canonical moment of ^ defined in Definition \2. 3[ and H^+i is the {m+l)th 
Verblunsky coefficient defined by the Szego recursions fl3.21l) and (I3.22I) . then 



(3.23) 



A, 



m+l — H^a+l 



Proof: Integrating the recursion fl3.22p we obtain 
and 

where we have used the orthogonality of the matrix polynomials stated in Theorem 13. II 

and the representations of the orthonormal polynomials 0^ and 0^ in Corollary 13. 2[ Observing 
Theorem 13.11 and the identity 



(3.24) 

yields 
(3.25) 



J —n 



m+l 



r-l/2/r ^L\-l/r ^O/flX pl/2 
-^m Vpi ^ ml R Vpj m/ R-'-^m 

T-l/2lnn 1-1 /r ^^,R\ pl/2 

-l^m l-'-ml \^p^^^mlR^m ■ 



The matrix polynomial '^^{z) has the representation 



m—l 



R^k 
k ^ 5 



A:=0 



where K^, . . . ,K^_^ denote p x p matrices and the leading coefficient is given by fl3.18p . 
Integrating with respect to dfi{9) gives 



k=0 
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and it follows from fl3.25p that 

m— 1 

(3.26) H^^, = L-'/'{I„ + ^ [L^y' z'^')nR-J". 

k=0 

Observing the definition of the canonical moments in (12.171) and the definition of the center (12.151) 
the assertion of the Theorem follows if the identity 

m—l 

(3.27) {i„z"^+' + Y,K^{L^y'^'''')R = r„+i-(ri...r^)T-ii(r_^...r_i)*. 

fc=0 

can be established. For this purpose we determine the matrices {k = 0, . . . ,m — 1) explicitly 
using the representation of the orthogonal matrix polynomials "^mi^) (13.61) . From this definition 
it follows that the element in the position of the matrix is obtained by deleting the 
(mp + j)th row and the {kp + i)th column in the determinant T^^{z) defined by (13.81) . that is 

Here again denotes the matrix obtained by deleting the {mp + j)th row and 

{kp + i)th column, which coincides with the entry in the position {kp + i, mp + j) of the adjoint 
of the matrix T^. Consequently, it follows that 

and the "vector" 



\T 







coincides with the right upper block of size mp x p of the matrix T^^. By standard result in 
linear algebra this block is given by 



m y 



~Tm-l(^~m ■ ■ ■ r_i)*_R, 

which yields 

m— 1 m—l 
k=0 k=0 

= -l^rnKFl . . . Tm)T^-i(r~m ■ ■ ■ r_l)*-R„^ 

Combining this result with the identity [L^^ ^ = Rm\Tm\~^ finally gives (13.271) . which completes 
the proof Theorem 13.31 □ 
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4 Geronimus relations for monic polynomials 



In this section we present a new proof of the Geronimus relations, which provide a representation 
of the canonical moments (or Verblunsky coefficients) of a symmetric matrix measure on the 
unit circle in terms of the coefficients in the recurrence relations of a sequence of orthogonal 
polynomials with respect to an associated matrix measure on th e interval [—1,1]. There exists 
seve ral alternative proofs of these relations in the literature [see lYakhlef and MarcellanI (120011 ) 
(120081 )] ■ but the one prese nted here explicitly uses t he theory of canonical 



and 



Damanik et al. 



moments of matrix measures as introduced in 



Dette and Studded (120021 ) . As a by-product we 



derive several interesting properties of the Verblunsky coefficients. 

To be precise let fic denote a symmetric (with respect to the point 0) matrix measure on the 
unit disc (i.e. fic is invariant with respect to the transformation 6 i-^ —0)- We associate to 
fic a corresponding matrix measure, say /i/, on the the interval [—1,1], which is defined by the 
property 

(4.1) [' f{x)dfij{x) = r /(cos {e))dfici0) 



for all integrable functions / defined on the interval [—1, 1]. Note that the relation Sz : dfiQ ^ dfij 
is called Szego mapping in the literature, where the matrix measure /i/ is usually defined on the 
interval [—2,2]. We will work wi th the i nterv al [— 1, 1] in this sectio n, because this interval is 
also used in the classic al papers of ISzegol (Il922[) a nd 
canonical moments by lDette and StuddenI (119971 ). 



Geronimus! (119461 ) and in the monograph on 



Note that the inverse of the Szego mapping (14. ip is characterized by the property 



(4.2) 



g{0)dfic{0) 



5f(arccos {x))dfij{x) 



-1 



where g denotes any integrable function on with g{6) = g{—9) for all 9 G [— 7r,7r). For a 
proof of the Geronimus relations we need several preparations. Our first results shows that the 
canonical moments (or Verblunsky coefficients) of a symmetric matr ix measure on t h e unit circle 
are real and symmetric matrices. The result was also proved by 



Damanik et al. 



mm . We 



provide here an alternative proof, because several steps in the proof are used later. 

Lemma 4.1 For any symmetric matrix measure nc on the unit circle the corresponding canonical 
moments are real and symmetric. 



Proof: By the symmetry of the matrix measure /i^ we have 



r.^e^^'df^cie) 



J^^e ^^^dfic{9) = r_fc which yields Tk = j^^cos{k9)diiciP). Consequently, the block Toeplitz 
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matrix associated with fic is given by 



/ Tn ... \ 



(4.3) 



T 



V 



r, 



/ 



and is symmetric. Because all entries of the matrix are real, the canonical moments are 
also real and it remains to establish the symmetry. 

For this purpose we denote by the p x p block in the position (fc, /) of the mp x mp— block 

matrix A. We will show at the end of this proof that 

[^m-l] (^k,l) ^ [^m-l] (m+l-fc,m+l-0 ' 

From this identity and the property = F^ we obtain 

m m 

(Fi, . . . , Tm)T^^i(Tm, . . . , Fi)* = ^ Ffc [r^ii] Tm+l^l = ^ Tm-k+1 [^m-l] 

k,l=l k,l=l 
m 

= ^ Tm-fc+l [^m-l] (fc^i) = (Xm, ■ ■ ■ , ri)T^ii(Fi, . . . , F^)*, 
k,l=l 

and by similar arguments 

(4.5) (Fi, . . . , Fm) T^i]^ (Fi, . . . , Fm) = (F^, . . . , Fi) T^i^^ (F^, . . . , Fi) . 



Observing the definition of the canonical moments A^+i it now follows that 

^m+i = — (Fm, . . . , Fi) T^i^ (Fm, . . . , Fi)*] ^ (F-m+i — (Fi, . . . , Fm) T^ii (Fm, . . . , 
X [Fo-(Fi,...,FjT,;ii (Fi,...,Fj*]-'/' 

= Am+l 

which proves the remaining assertion of Lemma 14.11 



ri)*)^ 



Proof of the identity 04.41) . The element in the position {i,j) of the matrix [T!^!^]^^^^ and 



{m+l—k,m+l—l) 



are given by 



and 



l/i \ {2m-l-k)p+i+j 



|T^_i|-^(-l) 



m— 1 



rri{im-l)p+j),({m-k)p+i) 
m—1 
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respectively, where T^'^i ''^^^^^'^^'^ ^jg^Qi^gg i^j^g matrix obtained from Tm-i by deleting the 

{m — l)p + j row and {m — k)p + i column (note that both expressions have the same sign). In the 
following discussion we denote by A^'^'^"^^ and A^^'^'^'^ the matrix obtained from A by deleting the 
2th column or the jth row, respectively. Then interchanging first columns and then rows yields 



T, 



m—1 



r 
r 



1-2 

l-l 



k-2 



\l-k\ 



r 



m—l 



r|«-fc+2| 



m—k+1 



-1)' 



-1) 



27 



Tm— 1 



J- m-Z+1 

^ m-l 
Tm-Z-l 



ro 

To 

r(j),{-) 

Tm— 1 



r ■ 

^ \l-k-l\ 

r(i),W 



m— fc 

Tk 



J- |«-fc-2| 

r(i),(-) 
^ \i-k-i\ 

^\i-k\ 



r 



m—k—l 
(■),« 



m—l 

r(j),(-) 



fc-2 



r 



\l~k~l\ 

r 



r(-),W 

^ \i-k~i\ 
r(i),» 



i \l-k+l\ ^ \l-k+2\ 



J- m— fe— 1 -L m—k m—k+1 

r{-),» r 

i ™ i m-k+l 



•m—k—l 



r|/-fc| 

^ \l-k-l\ 
'^\l-k-2\ 



m—k 



r{-),» 
r{i),W 
{■),« 



r 



r 



{■),» 

fc-i 



r|«-fc+2| 

r(i),{-) 

r|i-fc| 



rfc-2 



r;-2 



1 

Tm— 1 



>{■: 

L 

ri-2 



r(i),{-) 



{{m-l)p+j),{{m-k)p+i) 



m—l 



for some 7 G N, because the number of changed columns coincides with the number of changed 
rows. This implies (14.41) an completes the proof of Lemma 14.11 □ 



For the next step we need to define canonical moments of matrix mea sures on the interval [— 1, 1 ]. 
Because the main arguments here are very similar to the proceeding in lDette and StuddenI (120021 ). 
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who considered matrix measures on the interval [0, 1], we only state the main differences without 
proofs. To be precise, define for a matrix measure //,/ on the interval [—1,1] the moments Sk = 
Skil^i) = /^i x'^d/iiix) {k^O,l,.. .) and a vector Cn{iii) = {SoM, SnM) G {Wp)^+\ We 
consider the moment space 

(4.6) = {cnM I /// is a matrix measure on [-1, 1]} C (RPXf)"+i 

corresponding to the first n moments of matrix measures on the interval [—1, 1]. For a matrix 
measure /ij on the interval [—1,1] we define the block Hankel matrices Hj and I^j 

( So ... \ 



Sm ... S' 



2m j 



^ So — So 



H 



2m 



ii-2m+l 



H 



2m+l 



We introduce the notation 



Y Sm-l S„i+1 

I So + S^ . 

y Sm + 'S'm+l ■ 
/ ^0 - ^1 . 

Y Sm ~ Sm+l 



Sm—1 Sm+l ^ 



s 



2m-2 — S2m j 



Sm + Sm+l 



\ 



S2m + S2m+1 J 
Sm ~ Sm+l ^ 

S2m — S2m+1 J 



ihm ~ {Sm: ■ ■ ■ ■> S2m-l)'^ : ^2m — {Sm-l ~ Sm+l, ■ ■ ■ , S2m-3 ~ S2m-l)^ , 



h. 



■2m+l 



{Sm + Sm+l, ■ ■ ■ , S2m-1 + S2m)'^ , ^2m+l — {Sm ~ Sm+l, ■ ■ ■ , S2m-1 ~ S2m)'^ , 



and define S^ = Sq, S2 = Sq, 



S 



(4.7) 

and S^ = —So, 
(4.8) 



2m 



Stm+l 



■'2m 



^2m+l 



S2m-2 - ^2m-f^2m-2^2m {m>2), 



S2m — ^2m+l-f^2m-1^2m+l {m > 1), 



h^mH2m-2h.2m > 1), 

I^m+lS.2m-lh.2m+l " S2m {m > 1). 
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Note that the q uantities S:^ and S*" a re det ermined hj Sq, . . . , Sn-i- It can be shown by the same 
argument as in iDette and StuddenI (120021 ) that for {Sq, . . . , Sn-i) € Int(A^„) and any matrix 
measure fij on the interval [—1, 1] with moments satisfying Sj{ni) = Sj {j = 0, . . . ,n — 1), the 
moment of order n SnifJ-i) = j\x"'dfij(x) satisfies 



(4.9) 



With these preparations we can define the canonical moments of a matrix measure on the interval 
[— 1, 1] with moments 5*0, ... , Sn-i- 

Definition 4.2 Let fij denote a matrix measure on the interval [—1, 1] with moments Sk = 
SkifJ'i) = j\x^diJii{x) (^A; = 0, 1, . . .) and define 



(4.10) 



iV(/i,) = min |A; G N I (5o, . . . , 5^) G dM^^l^] 



For any n = 0, . . . , N{hj) — 1 the (symmetric) canonical moments of the matrix measure ni are 
defined by 



(4.11) 



-1/2 



where the quantities S^^^ and S^^^^ are given by ^^"-^ U-^ j respectively. 



Note that iDette and Studded (120021 ) use a non symmetric definition of canonical moments of 
matrix measures on the interval [0, 1], that is 



(4.12) 



Un+l — {Sn+l S^_^_i) {Sn+1 S^_^_i) . 



This non symmetric definition turns out to be more useful when working with monic orthogonal 
polynomials but in the present context the symmetric version has advantages. We are now in 
a position to prove the main result of this section, which relates the canonical moments of a 
symmetric matrix measure on the unit circle and the canonical moments of the associated matrix 
measure on the interval [—1, 1] by the Szego mapping. For this purpose recall the definition of 
the matrix ball in (12.141) and the defintion for the matrices L^, and (I2.12p . (12.131) 
and (12.151) . respectively. If the given measure /i^ on the unit circle is symmetric, then it follows 
from (145|) 



(4.13) 



The following result is the main step for the proof of the Geronimus relations. 
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Theorem 4.3 Let /ic denote a symmetric matrix measure on the unit circle and denote by 
fii = Sz{fic) the associated matrix measure on the interval [—1, 1] defined by the Szego mapping 
Iji4-1\ )- The canonical moments An and Un of the matrix measures fic O'^d /i/ satisfy 

An = 2Un-Ip; n = l,...,N{fic). 
Similarly, the non symmetric canonical moments Un defined in Iji4-^^ satisfy 

(4.14) 2Un-Ip = An; n = l,..., N{fxc), 
where the quantities An are given by 

(4.15) An = L-'_{'AnLi;\. 

Proof: We only prove the first part of the Theorem. The second part is shown by similar 
arguments. Assume that m < N{fic) and let Fq, Fi, . . . , denote moments of the matrix measure 
on the unit circle /ic. For j = 0,1,... we define Tj{x) = cosf? arccosx) as the ?'th (scalar) 



Chebychev polynomial of the first kind, then it follows from (14. 2 p and from iRivlinI (jl990l ) that 

/7r pi 
COS {j9)d^c{d) = / Tj{x)dfii{x) 

Li/2J 



where Si = j\x^diii{x) {I = 0,1,...) denote the moments of the associated matrix measure 
/i/ = Sz{nc) on the interval. Recall the definition of Sj^+i and S^j^^ in (14.71) and (14. Sp . then there 
exist matrix measures /i/ and /i7 on the interval [—1, 1] such that Sj = Sj{fif) {j = 0, . . . , m) 
and _^ ^ 

^++,=^^x-+M/x|(x) and S^^i = ^^x"+M/i7(x). 

We define 

(A T7\ r+ - 9"^Q+ _L \^ r ("^ + l)r(m + 1 - fc) _2fco 
(4-17) r^^, - 2 Sn,^, + 2. (-1) r(A: + l)F(m-2A: + 2)^ 

k=l 

r4 1Q^ r- - 9"^^- , (m+l)F(m + l-A;) ^^,2^ 

(4-18) r„+i - 2 ^ i-l) r(fc + l)F(m-2A; + 2)^ 

With the inverse Szego mapping we obtain the symmetric measures = {Sz)^^{fif) and fi^ = 
{Sz)^^{fij) on the unit circle and the representation ( ]4.16p yields that the measures and fi^ 
satisfy 

/n pn 
cos{{m + l)e)d^i+{0)=Ti+^ and / cos ((m + l)^)(i/Xc(^) = T'+i. 
-7T J —TV 



20 



Consequently, recalling the definition of the set Km in fl2.14p we have T^^^, F^^^ G and from 
the extremal property of the moments and 5*^^+1 we obtain that r^^^,r~^^ G dKm- By 

the definition of the set Km in (12.141) it therefore follows that the canonical moments A^+i 
Am+i corresponding to matrix measures /ij and /i^, respectively, are unitary. Moreover, Lemma 
I4.1[ implies that the matrices A^^^ and A^^^ are symmetric with real entries, which yields 

(^m+l) = i^m+l) = ^P- 

Consequently all eigenvalues of the matrices A^+i ^m+i given by —1 and 1. 
We now define the matrices 



(4.19) f++i = M„ + L^ and T^+i = - L 



mi 



which are obviously elements of the set Km because by (I4.13P we have Lm = Rm- Consequently, 
there exist matrix measures jl^ and jl'^ such that Tj(jl^) = Tj (j = 0, . . . , m) and 

Without loss of generality we assume that /ij and fl^ are symmetric with respect to the point 
[otherwise use |(/ij(6') + /ij(— 6*))] and we define fif = Sz{fi^) and ji^ = Sz{fl'^) as the associated 
measures on the interval [—1,1] with (m + l)th moments 5*^+1 and Sm+i, respectively. These 
matrices satisfy the identities 



f+ _ 9m Q+ I \^ / T\k {m + l)T{m + 1 — k) ^m~2kQ 
Im+i - ^ + 2^ ( 1) Y(k + l)T(m-2k + 2r 



k=l 

L(m+1)/2J 



_ om$- , V- f ,^k im+l)T{m + l-k) ^m^^kc 
- 2 bm+i+ 2^ ( r(A; + l)r(m-2A; + 2)^ 



k=l 



From the inequalities (14. 9 p it follows that Sm+i > 'S'j^+i and Sm+i > S~_^_i (note that 5'^_|_i and 
5*^^^ are moments of a matrix measure on the interval [—1, 1] with moments Sq, . . . , Sm)- On 
the other hand we have 



{ Q+ _ Q+ \ - -r+ _ r+ 



m+1 

Mm + Lm- (Mm + L]1^ A^^^L]!,^ 



— Vp ^m+l) 

> 0, 
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because the eigenvalues of the matrix Ip — Am+i are given by and 2. So we obtain 



while a similar argument shows 
Consequently, it follows that 



C- _ c- 



A. 



m+l 



m+1 



•-P ) 



A~ — —T ■ 



r+ 

m+l 



F" = r~ 

^ m+l ^ m+l ' 



and we obtain from the definitions of T^^^, T^_^^ in fl4.19p 

Mm = ^i^m+l ~^ ^m+l)^ ^rn = '^("'^m+l ~ '^m+l) 



The definition of the (m + l)th canonical moment Am+i of the matrix measure /i and fl4.17p - fl4.18l) 
now imply 



-1/2 



Am+1 = Lj^^'^{rm+1 — Mm)L^^'^ 

~ {ij^^m+l^'^m+l)^ rn+1 ~ -^(j^m+1 ^ m+lj^ i^'^^m+l ~ m+lj 

— {^m+1 ~ 'S'm+l) (2S'm+l — (S^+l + S^+l)) {S^n+1 ~ ^m+l) 

~ 2 — ^m+l) [pm+l — ^m+l) ('^'m+l ~ "^m+l) '~ 

= 2Urn+l — Ip, 

where the last equality is a consequence of the definition of canonical moments of matrix measures 

on the interval [—1, 1]. This proves the assertion of the theorem. □ 

Our final result gives the Geronir nus relatiori s for m onic orthogonal matrix polynomials, w hich 
generalize the results obtained by iGeronimud fll946l ) and iFaybusovich and GekhtmanI fjl999l ) for 
the scalar case. To be precise note that Corollary 3.2 together with (14.131) yield for the monic 
orthogonal polynomials and defined in (13.141) and (13.151) . respectively 



L iL 
Pm'Pm+l 

4>l 



r'l/2^L ,R R ^ J -1/2 

m m+l) Tm+lrm m+l m 

l^-l/2^R IL ^ j^-1/2 
m m' Vm m m 



Using these equations we obtain from fl3.2ip . fl3.22p and the second part of Theorem 14. 31 the Szego 
recursion for the monic orthogonal matrix polynomials with respect to a matrix measure on the 
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unit circle, that is 

)^ ( 



Consequently, the matrices A^+i defined by f l4.15p are the Verblunsky coefficients corresponding 
to the monic orthogonal polynomials and we obtain the following result. 

Theorem 4.4 Let fic denote a symmetric matrix measure on the unit circle and denote by 
fii = Sz{fic) the associated matrix measure on the interval [—1, 1] defined by the Szego mapping 
Iji4-1\ )- If Po> Pi! - ■ ■ be the monic polynomials orthogonal with respect to the matrix measure /ij 
satisfying the three term recurrence recursion 

(4.20) (1 + t)P„+l(t) = Pm+2it) + Rm+l{t)Cm+l + Pm{t)Bm, 

(Po{t) = Ip, -P-i(t) = Op), then the matrices B„i and Cm+i satisfy 

1 — —2 — 

where the quantities An are defined in fl4.15p . 



Proof: It follows analogously to iDette and StuddenI (120021 ) that the matrices Bm and C^+i are 
given by 

Bm = {S2m — S2m) {S2m+2 ~ S2'm+2)^ 
Cm+1 = {S2m+2 — S2m+2) ^{S2m+3 — S2m+3) + {^2171+1 " 'S'g^+i) "^(5'2m+2 " S2m+2)- 

and that the non symmetric canonical moments defined by fl4.12p satisfy 

"^Vn-lU n = {Sn~l S^—i) {Sn 'S'^ ), 

whenever n < N{fij), where Vn = Ip — Un- Consequently, the assertion follows by a direct 
application of the second part of Theorem 4.3. □ 
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